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I/"") i Dynamical chiral symmetry breaking is a nonperturbative phenomenon that may 

VP , be studied using QCD's gap equation. Model-independent results can be obtained 

C " ) ■ with a nonperturbative and symmetry preserving truncation. The gap equation 

yields the massive dressed-quark propagator, which has a spectral representation 
C _ ) when considered as a function of the current-quark mass. This enables an expli- 

, cation of the connection between the infrared limit of the QCD Dirac operator's 

C "3 ' spectrum and the quark condensate appearing in the operator product expansion. 



1. Introduction 

o ; 

3 . Dynamical chiral symmetry breaking (DCSB) is fundamental to the strong 

interaction spectrum. It is the generation through self-interactions of a 
momentum-dependent quark mass, M(p 2 ), in the chiral limit, that is large 
in the infrared but power-suppressed in the ultraviolet: 

m(/) lars = p2 ^ Hgg)°) (1) 



p2|ln 



A 2 

QCD 



The phenomenon is impossible in weakly interacting theories. In Eq. (1), 
7 m = 12/(33 — 2Nf) is the mass anomalous dimension, with Nf the num- 
ber of light-quark flavours, and (qq)° is the renormalisation-group-invariant 
vacuum quark condensate. 1 Dynamical chiral symmetry breaking is an es- 
sentially nonperturbative effect. 

It is natural to explore DCSB using QCD's gap equation and the evo- 
lution of the dressed-quark mass- function in Eq. (1) to a large and finite 
constituent-quark-like mass in the infrared, M(0) ~ 0.5 GeV, is a long- 
standing prediction of such Dyson-Schwinger equation (DSE) studies 2 that 
has recently been confirmed in simulations of quenched lattice-QCD. 3 The 
features and flaws of the approach are well known. For example, the DSEs 
are a keystone in proving renormalisability and provide a generating tool for 
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perturbation theory. However, they are also a collection of coupled integral 
equations from which a tractable problem is only obtained by truncation. 
A weak coupling expansion, which reproduces perturbation theory, is one 
systematic truncation procedure. However, that cannot be used to study 
nonpcrturbative phenomena. 

2. Truncating the DSEs 

Fortunately a systematic alternative is available. 4 Its leading order term fur- 
nishes a reliable description of vector and flavour nonsinglet pseudoscalar 
mesons and their interactions, 5,6 and the reasons for success in these chan- 
nels can be understood a priori. The scheme is worth illustrating. 
The renormalised gap equation 

Sip)' 1 =Z 2 (t 1 -p + m hm )+ Z, j\ 2 D^(p- q)^S(q)T a u (q,p) , (2) 

involves: D^Qt), the dressed-gluon propagator; T®(q;p), the dressed- 
quark-gluon vertex; m bm , the A-dependent current-quark bare mass that 
appears in the Lagrangian; :— j A cftq/ {2ir) i , a translationally-invariant 
regularisation of the integral, with A the regularisation mass-scale; and 
the quark-gluon- vertex and quark wave function renormalisation constants, 
Zi(( 2 ,A 2 ) and Z 2 {C > ' 2 1 h?) respectively, which depend on the renormalisa- 
tion point, the regularisation mass-scale and the gauge parameter. 

The kernel of Eq. (2) is formed from a product of the dressed-gluon 
propagator and dressed-quark-gluon vertex. It is to this product that a 
systematic truncation must be applied. However, one must also consider 
more than just the gap equation's kernel. Chiral symmetry is expressed via 
the axial-vector Ward-Takahashi identity: 

P l >T^{k-P) = S- 1 {k + )i lb +i lb S- 1 (k_), (3) 

k± — k ± P/2, wherein T^^k'jP) is the dressed axial-vector vertex. This 
three-point function satisfies an inhomogeneous Bethe-Salpeter equation: 

[IV(fc; P)] tu = Z 2 [ l5l ,] tu + [ A [S(q + )T 5 ,(q; P)S(q.)] sr K^(q, k; P) , (4) 

in which K(q,k;P) is the renormalised fully-amputated quark-antiquark 
scattering kernel, and the colour-, Dirac- and flavour-matrix structure in the 
equation is denoted by the indices r,s,t,u. The Ward-Takahashi identity, 
Eq. (3), means that the kernel in the gap equation and that in the Bethe- 
Salpeter equation are intimately related. Therefore a qualitatively reliable 
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understanding of chiral symmetry and its dynamical breaking can only be 
obtained using a truncation scheme that preserves this relation, and hence 
guarantees Eq. (3) without a fine-tuning of model-dependent parameters. 

The truncation scheme introduced in Rcf. [4] is a dressed-loop expan- 
sion of the dressed-quark-gluon vertices that appear in the half-amputated 
dressed-quark-antiquark scattering matrix: S 2 K , a renormalisation-group 
invariant. 5 All n-point functions involved thereafter in connecting two par- 
ticular quark-gluon vertices are fully dressed. The effect of this truncation 
in the gap equation, Eq. (2), is realised through the following representation 
of the quark-gluon vertex, r" = ^X a T^: 



Here T 3g is the dressed-three-gluon vertex and it is apparent that the lowest 
order contribution to each term written explicitly is 0(g 2 ). The ellipsis 
represents terms whose leading contribution is 0(g 4 ); viz., the crossed- 
box and two-rung dressed-gluon ladder diagrams, and also terms of higher 
leading-order. 

This expansion of S 2 K, with its implications for other n-point functions, 
yields an ordered truncation of the DSEs that guarantees, term-by-term, 
the preservation of vector and axial- vector Ward-Takahashi identities, a 
feature that has been exploited 7,8 to prove Goldstone's theorem and other 
exact results in QCD. It is readily seen that inserting Eq. (5) into Eq. (2) 
provides the rule by which the rainbow-ladder truncation can be systemat- 
ically improved. 

With this scheme, as with perturbation theory, it is impossible, in gen- 
eral, to obtain complete closed-form expressions for the kernels of the gap 
and Bethe-Salpeter equations. However, for the planar dressed-quark-gluon 
vertex depicted in Fig. 1, closed forms can be obtained and a number of 
significant features illustrated 5 when one uses the following model for the 
dressed-gluon line 9 



Z 1 r^(k, P )=^ + -^- [ g 2 D pa {p-l) lp S{l + k-p) 1 ,S{i) la 



N f 

+ -£ / g 2 D a , a (t)D T , T (e + k-p) 7T ,S(p-t)^rl%(£,-k,k-p) 




(5) 



g 2 D iiV (k)= 5^- 



kft k, 



(2tt) 4 G 2 5\k), 



(6) 



where G sets the model's mass-scale. This model has many positive features 
and, furthermore, its particular momentum-dependence works to advantage 
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Figure 1. Integral equation for a planar dressed-quark-gluon vertex obtained by neglect- 
ing contributions associated with explicit gluon self-interactions. Solid circles indicate 
fully dressed propagators. The vertices are not dressed. (Adapted from Ref. [5].) 

in reducing integral equations to character-preserving algebraic equations. 
Naturally, there is a drawback: the simple momentum dependence also 
leads to some mo del- dependent artefacts, but they are easily identified and 
hence not cause for concern. 

It is a general result 5 that with any vertex whose diagrammatic content 
is known explicitly; e.g., Fig. 1, it is possible to construct a unique Bcthc- 
Salpeter kernel which ensures the Ward-Takahashi identities are fulfilled: 
that kernel is necessarily nonplanar. This becomes transparent with the 
model in Eq. (6), using which the gap equation obtained with the vertex 
depicted in Fig. 1 reduces to an algebraic equation, irrespective of the 
number of dressed-gluon rungs that are retained, and the same is true of 
the Bcthc-Salpeter equations in every channel: pseudoscalar, vector, etc. 

Results for the tt and p are illustrated in Table 1. It is apparent that, 
irrespective of the order of the truncation; i.e., the value of n, the number 
of dressed gluon rungs in the quark-gluon vertex, the pion is massless in 
the chiral limit. (NB. This pion is composed of heavy dressed-quarks, as 
is evident in the calculated scale of the dynamically generated dressed- 
quark mass function: M(Q) « 0.5 GeV.) The masslessness of the tt is a 
model-independent consequence of consistency between the Bethe-Salpeter 
kernel and the kernel in the gap equation. Furthermore, the bulk of the 

Table 1. Calculated tt and p masses, in GeV, quoted with G = 0.48 GeV. n is the 
number of dressed-gluon rungs retained in the planar vertex, see Fig. 1, and hence the 
order of the vertex-consistent Bethe-Salpeter kernel: the rapid convergence of the kernel 
is evident from the tabulated results. (Adapted from Ref. [5].) 
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0.695 


0.762 


0.770 


0.770 



February 8, 2008 23:10 WSPC/Trim Size: 9in x 6in for Proceedings 



cdrqc5 



5 

p-ir mass splitting is present in the chiral limit and with the simplest (n = 
0; i.e., rainbow-ladder) kernel, which shows that this mass difference is 
driven by the DCSB mechanism: it is not the result of a finely adjusted 
hypcrfine interaction. Finally, the quantitative effect of improving on the 
rainbow-ladder truncation; i.e., including more dressed-gluon rungs in the 
gap equation's kernel and consistently improving the kernel in the Bethe- 
Salpeter equation, is a 10% correction to the vector meson mass. Simply 
including the first correction (n = 1; i.e., retaining the first two diagrams 
in Fig. 1) gives a vector meson mass that differs from the fully resummcd 
result by < 1%. 

While these results were obtained with a rudimentary interaction model, 
the procedure is completely general. However, the algebraic simplicity of 
the analysis is naturally peculiar to the model. With a more realistic in- 
teraction, the gap and vertex equations yield a system of twelve coupled 
integral equations. The Bethc-Salpeter kernel for any given channel then 
follows as the solution of a determined integral equation. 

By identifying the rainbow-ladder truncation as the lowest order in a 
systematic scheme the procedure also provides a means of anticipating the 
channels in which that truncation must fail. The scalar mesons are an ex- 
ample. Parametrisations of the rainbow-ladder truncation, fitted to ir, p ob- 
servables, yield scalar mesons masses that are too large. 10 That was thought 
to be a problem. However, we now know this had to happen because can- 
cellations that occur between higher order terms in the pseudoscalar and 
vector channels, thereby reducing the magnitude of corrections, do not 
occur in the scalar channel, wherein the full kernel contains additional 
attraction. 11 Therefore an interaction model employed in rainbow-ladder 
truncation which simultaneously provides a good description of scalar and 
pseudoscalar mesons must contain spurious degrees of freedom. Quanti- 
tative studies of the effect of the higher-order terms have begun 11 ' 12 and 
a straightforward understanding of scalar mesons is not impossible. Simi- 
larly, it may be that the only thing remarkable about mesons with "exotic" 
quantum numbers is that their accurate description requires only an equally 
carefully considered application of this systematic approach. 13 

3. Concerning the quark condensate 

While Eq. (1) is expressed in Landau gauge, the so-called OPE condensate, 
(qq) , is gauge parameter independent. In the chiral limit this conden- 
sate sets the scale of the mass function in the ultraviolet and thus plays 
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a role analogous to that of the renormalisation-group-invariant current- 
quark mass in the massive theory. A determination of the OPE conden- 
sate directly from lattice-QCD simulations must await an accurate chiral 
extrapolation 14 but DSE models tuned to reproduce modern lattice data 
give 15 |(gg) |~A3 CD . 

Another view of DCSB is obtained by considering the eigenvalues and 
eigenfunctions of the anti-Hermitian massless Euclidean Dirac operator: 

7 • Du n (x) = iX n u n (x) . (7) 

The eigenfunctions form a complete set, and except for zero modes they oc- 
cur in pairs: {u n (x) , j$u n (x)} , with eigenvalues of opposite sign. It follows 
that in an external gauge field, A, one can write the quark Green function 

S(x, y; A) = (q(x)q(y))^ = £ , (8) 

~r III 

n 

where m is the current-quark mass. Assuming, e.g., a lattice regularisation, 

i//*«K*), (l »s = -f £ (•) 

where V is the lattice volume. a One may now define a quark condensate: 

<0|gg|0r := Jim 1 / d*x { (q(x)q(x))% ) , (10) 

wherein the r.h.s. expresses an average over gauge field configurations. 
When V — > oo the operator spectrum becomes dense and Eqs. (9), (10) 
become 

-(0\qq\0) m = 2m f°° d\ , (11) 

with p(X) a spectral density. This equation expresses an assumption that 
in QCD the full two-point massive-quark Schwinger function, when viewed 
as a function of the current-quark mass, has a spectral representation. 
It follows formally from Eq. (11) that 

-<0|gg|0>°:= lim 2m [°° dX = n p(0) . (12) 

This is the so-called Banks-Casher relation, 17 which has long been advo- 
cated as a means by which a quark condensate may be measured in lattice- 
QCD simulations. 18 



a In deriving Eq. (9), zero modes have been neglected, which is justified under broad 
conditions. 16 
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The explication 19 of a correspondence between the condensate in Eq. 
(1) and that in Eq. (12) requires care because Eq. (11) is meaningless as 
written: dimensional counting reveals the r.h.s. has mass-dimension three 
and since A will at some point be greater than any relevant internal scale, 
the integral must diverge as A 2 , where A is the regularising mass-scale. 

The Schwinger function 

a(m) :=N c ix D f S m (p) , (13) 

which is the trace of the unrenormalised massive dressed-quark propagator 
evaluated at a fixed value of the regularisation scale, A, can be identified 
with the l.h.s. of Eq. (11). The essence of the Banks-Casher relation is that 
<T (to) have a spectral representation: 

t { rn),= 2rnj\x^, (14) 

where to = TO bm (A), and Eq. (14) entails 

~p{\) = ±- km [a(i\ + r?) - a(iX - rfi ] . (15) 

The content of this sequence of equations is elucidated by inserting the 
free quark propagator in Eq. (13). The integral thus obtained is readily 
evaluated using dimensional regularisation: 



O-free(m) = ^ m 3 



TO 1 

In —7- H V 7 - m 4n 

C 2 £ 



(16) 



With Eq. (15) the regularisation dependent terms cancel and one obtains 

In perturbation theory it can be shown that every contribution to p(A) is 
proportional to A 3 and hence 

(0\qq\0)° cx p{\ = 0) = 0. (18) 

A nonzero value of p(0) is an essentially nonperturbative effect. 

A precise analysis requires that attention be paid to renormalisation. 
Therefore consider the gauge-parameter-indcpcndent trace of the renor- 
malised massive quark propagator: 

<7(m; C) = Z 4 (C, A) A^ c tr D f S m (p; Q , 
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where the argument remains m = m bm (A). The renormalisation con- 
stant Z4 vanishes logarithmically with increasing A and hence one still 
has ct(to) ~ A 2 m bm (A). However, since 7 

,A 

-(qq)° ( = lira Z 4 (C,A)N c tv D S (p;C), (19) 

A^oo J p 

it is clear that for any finite but large value of A and tolerance S, it is always 
possible to find m bm (A) such that 

a(m; Q + {qq}° ( < 5 , Vm bm < m bm . (20) 

This is true in QCD and can be illustrated using the model of Ref. [1]. 
Figure 2 displays a(m,(), evaluated using a hard cutoff, A, on the inte- 
gral in Eq. (19), calculated with the massive dressed-quark propagators 
obtained by solving the gap equation. Since Eq. (20) specifies the domain 
on which the value of a(m; £) is determined by nonperturbative effects, one 
anticipates 

m bm (A) « -(qq)°/A 2 ~ 10~ 9 (21) 

for A = 2.0TeV in QCD where |(gg)°| ~ A^ CD . This estimate is confirmed 
in Fig. 2. 

The dotted line in Fig. 2 is 

ct(to,C) = 

-(qq)° c - arctan — + Z 4 ((, A)^m [A 2 - m 2 In (l + A 2 /m 2 )l , (22) 

and the difference between this and the curve is of order a (A) to (A) A 2 
because the DSE model incorporates QCD's one-loop behaviour. The fig- 
ure also displays <r(m, £) obtained in the absence of confinement, in which 
case 20 (qq)° = 0, as is apparent. 

It is now plain that in QCD a(m, () is a monotonically increasing 
convex-up function with a regular chiral limit and consequently it has a 
spectral representation: 

This lays the vital plank in a veracious connection between the condensates 
in Eqs. (1) and (12). On the domain specified by Eq. (21), the behaviour 
of <r(m,() in Eq. (19) is given by Eq. (22), which yields, via Eq. (15), 

7rp(A)^-(^ + Z 4 (C,A)^A 3 + ... , (24) 
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m bm (A) (GeV) 

Figure 2. Circles/solid-line: o-(rre) 1 '' 3 in Eq. (19) as a function of the current-quark 
bare-mass, evaluated using the drcsscd-quark propagator obtained in the model of Ref. 
[1]; dashed line: the model's value of (— {qq)® =1 q c v) = (0-24 GeV) 3 ; and dotted line: Eq. 
(22). Diamonds: cr(m) 1 / 3 evaluated in a non-confining version of the model; dot-dashed 
line: Eq. (22) with {gq)° = 0. (Adapted from Ref. [19].) 



where the ellipsis denotes contributions from the higher-order terms implicit 
in Eq. (22). 

Reference [21] reports the spectral density of the staggered Dirac oper- 
ator in quenched SU(3) gauge theory calculated on a,V — 4 4 -lattice, in the 
vicinity of the deconfining phase transition at (3 > 5.6. A comparison 19 
with these results shows that while the effect of finite lattice volume is ev- 
ident early (for Aa > 0.1, a is the lattice spacing) the behaviour at small 
Xa is qualitatively in agreement with Eq. (24) and Fig. 2: a nonzero OPE 
condensate dominates the Dirac spectrum in the confined domain; and it 
vanishes in the deconfined domain, whereupon p(0) = and the perturba- 
tive evolution, Eq. (17), is manifest. To be more quantitative, for a lattice 
coupling (3 — 5.4, p(0)a 3 s=s 70, so that 

np(0)a 3 /V « (0.95) 3 . (25) 

The lattice spacing was not determined in Ref. [21] but one can neverthe- 
less assess the scale of Eq. (25) by supposing a <~ 0.3 fm ~ 0.3/Aqcd, a 
value not unreasonable for small lattices and small (3, wherewith the r.h.s. is 
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~ (3Aqcd) 3 - This is too large but not unreasonable given the parameters of 
the simulation, its errors and the systematic uncertainties in our estimate. 
Fitting the lattice data at /3 = 5.8, one finds p(X) oc A 3 on A < 0.1 but with 
a proportionality constant larger than that anticipated from perturbation 
theory; viz. Eq. (24). Some mismatch was to be expected, however, because 
at P — 5.8 one has only just entered the deconfined domain and near the 
transition boundary some nonperturbative effects are still relevant. It is a 
modern challenge to determine those gauge couplings and lattice parame- 
ters for which the data are quantitatively consistent with Eq. (24). 

In practice, there are three main parameters in a simulation of lattice- 
QCD: the lattice volume, characterised by a length L; the lattice spacing, 
a; and the current-quark mass, to. So long as the lattice size is large com- 
pared with the current-quark's Compton wavelength; viz., L >• 1/m, then 
dynamical chiral symmetry breaking can be expressed in the simulation. 
Supposing that to be the case then, as outlined above, so long as the lattice 
spacing is small compared with the current-quark's Compton wavelength 

7T p(0) w - (qq)° 1/a ,o«l/m«L, (26) 

where the (qq)y a is the scale-dependent OPE condensate [£ = A = 1/a 
in Eq. (24)]. This completes an explication of the connection between the 
condensates. Note, however, that the continuum analysis indicates that one 
requires am < (oAqcd) 3 if p(A = 0) is to provide a veracious estimate of 
the OPE condensate. The residue at the lowest-mass pole in the flavour- 
nonsinglet pseudoscalar vacuum polarisation provides a measure of the OPE 
condensate that is accurate for larger current-quark masses. 19 
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